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, Various features of the mass yields in heavy ion collisions are studied. The mass 
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yields are discussed in terms of iterative one dimensional discrete maps. These maps 
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are shown to produce orbits for a monomer or for a nucleus which generate the mass 

yields and the distribution of cluster sizes. Simple Malthusian dynamics and nan- 
's , 

linear Verhulst dynamics are used to illustrate the approach. Nuclear cobwebbing, 
attr actors of the dynamics, and Lyapanov exponents are discussed for the mass distri- 
bution. The self-similar property of the Malthusian orbit offers a new variable for the 
study of scale invariance using power moments of the mass distribution. Correlation 
lengths, exponents and dimensions associated with scaling relations are developed. 
Fourier transforms of the mass distribution are used to obtain power spectra which 
are investigated for a 1// behavior. 
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I. INTRODUCTION 



The purpose of this paper is to illustrate how simple pictures developed from soluble 
models can lead to useful methods in understanding the mass yields in heavy ion 

collisions. First, we will attempt to show how iterative one dimensional maps j|,[5]|| can 
characterize ensemble averaged behaviors of these yields. The language of one dimensional 
maps contains such ideas as cobwebbing, orbits and trajectories, attractors, basins of attrac- 
tion, Lyapanov exponents, non-linear dynamics and chaos P,BLffl,^|,^|^S,0]. The approach 
taken here will also suggest new ways of portraying the experimental data which go beyond 
any specific model discussed here and which may give further insights into the dynamics of 
the collision. One main result of this paper is that the mass yield can be generated by the 
orbit of a monomer or of a nucleus (as the initial point) in its approach to an attractor. 
Simple Malthusian dynamics with its associated geometric growth will be shown to charac- 
terize some features of the data considered in this paper. This geometric Malthusian model 
will be shown to be an approximate representation of the soluble models considered in Refs. 
[1-3]. A non-linear Verhulst map will also be considered to investigate signals for non-linear 
dynamics and limitations to unbounded geometric growth. 

In Refs. [1-3] we have developed an exactly soluble model of fragmentation and obtained 
an evolution equation |2| which can be reduced to a Fokker-Plank equation for a mean 
cluster distribution. The evolution equation relates the change of the mean cluster number 
of a given size to the mean cluster number of all the other sizes when the parameter x 
(introduced in Ref. [1] and discussed below) or equivalently the temperature changes [||. 
The solution of this evolution equation at a given x value determines the mean distribution 
at the corresponding temperature. We have also shown |||| that our x model has a quite 
similar structure to a Markov population process [0] . These properties of our fragmentation 
model enable us to relate this model with one dimensional maps for the mean behavior in a 
fragmentation process. The discrete cluster size in a nuclear fragmentation corresponds to 
the iterative index of a discrete map. 
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Then, secondly, we will discuss the scaling and self-similar behavior of the power moments 
|| that can be obtained from the mass yield produced in a nuclear collision. Since a simple 
linear geometric map generates orbits having self-similar properties, power moments of the 
orbit points also will have a scaling power law behavior. The self similar property of the 
geometric or Malthusian model will be shown to offer a new variable for the study of scale 
invariance and universality in the power moments of the mass yield. Moreover, some analytic 
results for the scaling properties of the soluble model considered in Refs. UHJ^H wn l De given 
and compared to the simpler geometric model. The scaling properties of the power moments 
with bin size have associated correlation lengths and exponents related to the dimension of 
the system as discussed in Ref. [8] and further developed here. 

Thirdly, this paper will investigate the Fourier spectrum of the fragmentation distribu- 
tion. The importance of the Fourier spectrum is the appearance of a power law behavior 
with frequency /, such as 1/ f@ behavior, with (3 an exponent which, in a description of noise 
spectra, characterizes the type of noise ||. For example /3=0 is white noise, and (3=1, 1/ / 
noise. Finally, we will investigate some experimental data of Ref. []nj . 

II. GEOMETRICAL ITERATIVE MAPS AND MEAN CLUSTER 

DISTRIBUTIONS 

The interrelationships between one dimensional maps and mass yields developed out of 
a previous model of fragmentation and its generalizations Specifically, the model 

of Refs. |]|||§ contains a simple expression for the ensemble averaged mean distribution of 
cluster sizes Y^k, x) =< n k > given by 

_ x A\ T(x + A - k) 
<nk> ~k(A-k)\ T(x + A) ' () 

where the T's are gamma functions. The A is the total number of nucleons, k the size of 
the cluster, and the a; is a tuning parameter which can vary from to infinity depending on 
the temperature. The parameter x contains the dynamical quantities associated with the 
fragmentation. Thermodynamic arguments |],|2|,|3j give 
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V 

x = — exp 

v 

The V is the freeze out volume, T is the equilibriation temperature, and v Q is the quantum 
volume v = /i 3 /(27rm p fcsT) 3//2 with m p the mass of a nucleon. The a v is the coefficient in a 
simplified equation for the binding energy Eb of a cluster of k nucleons: Es{k) — a v (k — 1). 
The e Q is the level density parameter related to the spacing of excited levels in a cluster and T D 
is a cut-off temperature for internal excitations. Further details for this relation can be found 
in Refs. |l]J^J^]. Here, we mention the following. Small x corresponds to a low temperature 
and/or small interaction volume and large x corresponds to a high temperature and/or large 
interaction volume. Also for small x, the < >~ x and thus the < n& > depend on the 
separation energy a v in the Boltzmann factor exp(— a v /kBT), For large x, x » A, the 
< rik >~ Ai^A/xY' 1 and thus < >~ exp[a v (k — 1)//cbT] = exp[EB(k)/kBT). In this 
limit < rik > depends on the Boltzmann factor in total binding energy. The < rife > given 
by Eq.(|l|) thus contains both the evaporation region dominated by separation energies and 
the multifragmentation region dominated by binding energy factors p). As a final remark 
on the physical significance of x, we note that the exponent —ksTTo/ 'e (T + T ), which 
also appears as a Boltzmann factor, is the temperature dependent part of the Helmholtz 
free energy F per particle at low temperatures divided by ksT 0. In particular, F = 
(3iVei?/5)(l — 5ir 2 (k bT / e f) 2 / 12) where €f is the Fermi energy which is related to the level 
spacing e D via e Q = Aep/n 2 . Note, for low T, T Q /(T + T D ) = 1 and the dominant term in the 
factor exp[(— a v /kBT) — kT/e ] is the separation energy part, ie., (—a v /ksT). For large T, 
the factor exp[— ksTTo/ t a {T + T )] in x gives rise to the exponential behavior of the internal 
partition function of a cluster evaluated in a Fermi gas approximation. Thus, the factor 
exp[— kTT /e (T + T Q )] evolves from an evaporation factor into a Boltzmann enhancement 
factor for internal excitations. This behavior is similar to the evolution of the separation 
energy into a factor for the binding energy of a cluster from the exp(— a^/fc^T) part of x. 
In the rest of this paper, we will keep unspecified tuning parameter. The above 

remarks show that this tuning parameter contains the underlying physical quantities: V, T, 
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e Q and a v . 

The simple x-model of Eq.([I]) fits the smooth behavior of the fragment distribution in 



various proton-nucleus and nucleus- nucleus collisions [O]. However this model does not 
describe all the detailed behavior of these distributions such as fission peaks and shoul- 
ders [II|,r2] • O ne wa Y °f overcoming this limitation is using a size k dependent parameter x k 



instead of a single parameter x. The size dependence of Xk may arise from a more realistic 
treatment of the binding energy by including surface and Coulomb energies, and shell effects 
which were not included in obtaining Eq. ([]]). However, when this is done, we loose the exact 
solubability of the model. Another way of overcoming this limitation is to seek an alternate 
methods of approach for fragmentation process. Here, we will show that the result of Eq. ([I]) 
has an interpretation in terms of iterative one dimensional maps and use this interrelation 
as a starting point for another type of description of the complex dynamics of the fragmen- 
tation process. Iterative maps have been used to model the behavior of complex systems in 
many other areas ||. 

The x model of Eq.([T]) has a recursive or iterative property. Namely, the total average 
masses in clusters of size k, which is given by m k = k < n k >, can be related to the total 
masses of clusters of size k + 1 by 

m k+1 = r A (k,x)m k . (3) 

Here we dropped the broken bracket representing an ensemble average of m k for simplicity 
in notation. The recursive coefficient r A (k, x) is simply 

^■'> V1+*-i) BeMM ' (4) 

The initial point in the x-model is m\ = xA/(A + x — 1). The mass distribution then follows 
by successive iterations, 

mi — > m 2 — > m 3 — > • • • — > m k — > • • • — > m A , (5) 

and satisfies the constraint A = J2k m k- 



We can also consider the backward version of the forward invertible mapping given by 
Eq.©; 

m A _ k = r A {k,x)m A -k+i, 
r A (k,x) = l/r A (A - k,x). (6) 

Since Eq.@ is a linear invertible map, this backward map is the inverse map of the forward 
map and gives the sequence of 

m A -> m A _! y m k -> > m v (7) 

The sequence of Eq. (|5]) gives the orbit of m\ under iteration and the recursive relationship 
generates pairs (mk, mk+i) constructed by the recursive coefficient. This way of representing 
the results of the above model is similar to that used in studies of population dynamics and, 
in particular, to models in which the population evolves in discrete steps. The cluster size 
k serves as the discrete index and corresponds to the successive generations while the mass 
distribution rrik is the analogue of the population size in the k-th generation. The k runs 
from 1 to A — 1 and the orbit should be eventually normalized so that J2k=i m k — A. 

In a previous paper [i"3] another connection of the x-model of fragmentation with a 



model in population genetics due to Ewens |14| was explored. There the index k represented 



the number of times of a particular gene appears and n k the number of different types of 
genes each appearing k times in a sample of size A. The parameter x is connected with the 
mutation rate and effective population size in population genetics and evolutionary biology. 

Using the analogy with population dynamics, the data on the mass distribution can be 
plotted as points for each pair (m^, m^+i) for each k = 1, A — 1 with m^+i along the 
vertical axis and m^ along the horizontal axis. A smooth curve through the data gives the 
recursion curve for that data set. For, example, the simple Malthus model in population 
dynamics, which is 

""ifc+i = rmk = r k m\ = r fc+1 mo, (8) 
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gives a straight line for r 7^ 1 in such a plot. Here m = m-i/r and the constraint J2k=i m k = 
A determines mi = A(l — r)/(l — r A ). The parameter r can be determined by the data 
set by a least square fit and r represents the slope of a straight line fit. The parameter r 
can also be determined by the slope A = logr of a straight line fit in a logm^ vs k plot. In 
general this curve will be more complex than a straight line. It can depend on k as in Eq. (f|) 
for the x-model or can depend on rrik itself as in a non-linear map. The orbit of m\ upon 
successive iteration is determined by rrik+i = r k m,\ in the Malthus model. 

It should be noted that when the k dependence of r A (k, x) is averaged over then the frag- 
mentation model of Eq.(|3|) reduces to the simple Malthus model. For example, a geometric 
mean f can be defined by 

r = e x = [r A (l, x)r A (2, x) ■ ■ ■ r A {A - 1, (9) 
= exp ([X A (l,x) + X A {2, x) + ■ • • + \ A (A - 1, x)]/(A - 1)) . 

The A is the arithmetic mean of the Xyi{k,x) = lnr A (k,x). In the Malthus model or in 
the fragmentation model of Eq.(|3]) when r > 1 or x < 1, rrik increases geometrically with 
increasing k and when r < 1 or x > 1, decreases geometrically with increasing k. As 
we shall see, the mass distribution of cluster sizes in the back direction, i.e., k near A, is an 
important quantity. 

Iterative models also lend themselves to a geometrical iteration known as cobwebbing 
which involves the generation of the pairs (rrik, rrik+i) by graphical means. Starting with 
mi on the x-axis, is evaluated as the corresponding |/-axis value on the curve, y = rx, 
or m 2 = rm\. Then m 2 is reflected onto the x-axis through the 45° line, y = x. Next m 3 
is generated by repeating the procedure with m 2 on the x-axis and so on. The cobwebbing 
gives a useful visual way of understanding the trajectory or orbit and gives a large amount 
of information by pure geometrical means and simple graphical analysis. 

The behavior of rrik as k becomes large (k is limited by A in a nuclear fragmentation) gives 
the attractors of the orbits. In Malthusian dynamics, all orbits of m will go asymptotically 
to the +00 attractor or to the attractor if r > 1 or r < 1 respectively. The point m = is 
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an attracting fixed point for r < 1 and a repelling fixed point for r > 1. At r = 1 and from 
the result of Eq.(^) with x — 1, the orbit of m stays fixed for all values of m. For this case, 
rrik+i = m k forms a fixed line of the map. 

Figs. 1(a) - 1(d) illustrate an application of the above ideas to the experimental data in 
nuclear fragmentation. From the fits shown, we see that the fragmentation model of Eq. ([!]) 
with x = 0.3 (solid line) and its simple Malthusian version (dashed line) give a reasonably 
good fit to the data. 

These data show that the number of clusters decreases as the size k increases for small 
k, then slowly increases until it reaches a somewhat flat region (also called a "shoulder") at 
large k, then finally rises sharply again to k = A. The linear map of Eq.(H) fits the data 
better than the x-mode\ of Eqs.([|) and ([3D except for the p + Au data of Fig. 1(c). Allowing a 
simple k dependence in the coefficient r as r(k) = a + bk with a = 0.9833 and b = 0.000333, 
the linear Malthus map can fit the whole region of the p + Au data (dash-dotted line in 
Fig. 1(c)), (using a k dependent r is similar to mixing various r values as we will consider 
in Section VI for the Au emulsion data [|10|| ). This simple k dependent linear map has a 
geometric mean f of 1.016 and a minimum n k at k ~ 85. However both the x-model and 
the linear Malthus map do not produce a levelling off at a large mass region which may be 
related to a non-linear map as we will discuss in the next section. 



III. NONLINEAR MAP 

Non-linear models were introduced to limit the geometric growth of the Malthus model 
(attractor at oo). The classic examples are the Verhulst models; two examples are 

Rv m k , 1fV > 

and the logistic map, 

m k+1 = R v m k (K v - m k ) (11) 
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with the Ky — m k = if m k > Ky. The result of Eq.flllD can be rewritten as P k+ i = 
aP k (l — Pic)- The Ky is the carrying capacity of the system. The logistic map Eq.flTTD has 
two fixed points, m = and m = Ky — 1/Ry. Two fixed points of the Verhulst map of 
Eq.fllUD are m = and m = m s = Ry — Ky. These fixed points of Eqs.fllQ|) and (|11|) are all 
finite in contrast to the Malthus map Eq.(|S]). Since nuclei around k = 230 spontaneously 
fission, we expect a limit on the iteration number k, k max rs 230, when we consider nuclear 
fragmentation. It should also be noted that in the x-model the iteration is truncated after 
k = A — 1 since r^(A; = A, x) — 0. 

The Verhulst model given by Eq. (pTOj ) can be solved analytically because the inverse mass 
L k = obeys a linear mapping 

L k+ i = — + —^L k . (12) 

riy riy 

Due to this linearity, the Verhulst map has a unique inverse map and thus both the forward 
map and the backward map generate the same orbit once the end points (k = 1 and k — A) 
are the same. 

The resulting distribution of cluster sizes arising from the orbit of mi = n\ is given by 

n k +i = (t—t) n fc^^ k — ni ( 13 ) 

\k + 1/ (1 — Grjnx + a k m s 

where a = K v /R v and m s = Ry — Ky. The appearance of a k in Eq. flTBD gives a geometric 
behavior to n k+1 . When Ky > R v , L k — ► oo and m k — > geometrically. When Ky < R v , 
L k — > l/m s and m k m s , the saturating value (attractor) of the Verhulst mapping of 
Eq. (|10"D . Note n k falls because of the extra k in n k = m k /k. When Ry = Ky, then a = 1 
and n k+ i no longer has a geometric part but rather falls as a power law: 

Hk+1 = (kh) l + nlk/Ry ^ (M) 

For large k such that n\k/Ry » 1 the n k+ i ~ Ry/{k + l)k. 

The n k has the interesting property of having a bifurcation point when a = Ky / Ry = 
< 1. This behavior is not present in m k which increases uniformly with k for a < 1. 
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The orbit of monotonically approaches one of the two fixed points (0 or m s ) depending 
on the value of a. The orbit of can have both minimum and maximum points depending 
on the value of a and m s jm\. The bifurcation point arises from the extra k in going from 
mk —>■ nk = TThk/k. Specifically when m s jm\ > e 2_< » + 1, the nk has a minimum at k m i n 
and given by the solution to the transcendental equation y\ = y 2 with 

yi = ((k + 1) — 1 and y-i = e^ k /(m s /mi — 1). The value k m i n at which is a minimum is 
bounded from below by — 1 as m s jra\ — > oo. The k\/2 point is defined as mk 1/2 +i = m s /2 
and is given by k\/2 = [ln(m s /mi — 1)]/C At this point ku2, which is also the inflection 
point of rrik, the slope dyi/dk is equal to the slope dy2/dk. Thus k m i n < fci/2 < k max with 
kmin = ki/2 — k max when 

m s /m\ = e 2 ' c + 1. (15) 

Thus the condition Eq.(^) determines the bifurcation surface, i.e., the two roots k min and 
k max merge to k\/2 = 2/C, — 1, in the parameter space of (£, m s , and mi). 

Fig. 1(d) also illustrates a fit using the non-linear mapping of Eq.fllOl) to the p + Ag 
data where the maximum point k max is placed near the endpoint k = A to fit the shoulder 
(dash-dotted curve). Comparing the result of the x-model of Eq. ([[]), the Malthus model, 
and the Verhulst map near the endpoint with k ~ A, we see some differences. The result of 
Eq.(|I]) for x < 1 approaches the endpoint with an upward cusp, the Malthus model comes 
in as a straight line with positive slope for r > 1, and the Verhulst model can approach the 
endpoint with zero slope, positive slope, or negative slope depending on the position of k max 
with respect to the endpoint. The data show a flat region or "shoulder" near the endpoint 
and then finally a cusp rise. 

IV. SELF-SIMILARITY, SCALING BEHAVIOR AND UNIVERSALITY IN THE 

POWER MOMENTS OF MASS DISTRIBUTION 

In this section, we will study the behavior of the simple Malthus model of Eq.(^j) in 
greater detail. The trajectory behavior of a Malthus map can also be seen by looking at a 
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hypothetical two dimensional spiral 

m (r, 9) = m 1 r 9/2lT e i6 = x + iy, (16) 

where 9 can vary from to oo. The intersection of the continuous spiral with the positive 
real axis, i.e., m(r,9) = Xk at the values of 9 = 2-nih — 1), gives the discrete set of points 
x\ = m>i at 9 = 0, X2 = f«2 = rm\ at 9 = 2tt, and so on. The index k, representing the 
k-th intersection, is a labeling of clusters by size. The resulting points become the discrete 
"time" (t = k) version of a continuous system (t = 9/2n) in two dimensions. For r > 1 or 
x < 1, the system spirals outward to a fixed point at infinity and for r < 1 or x > 1, the 
system spirals inward to a fixed point at zero. For r = 1, the spiral reduces to a circle with 
a radius m\. In a fragmentation process, x = 1 corresponds to the point where the mass 
A is uniformly distributed over the cluster index k. The quantities A^(A;,z) which appear 
in Eq.(^j) can be considered as the Lyapanov exponents for each loop k and A = lnf the 
arithmetic mean Lyapanov exponent of the truncated spiral. 

We note that the geometric growth in m,k+i of Eq.(^) or m(r,9) of Eq.(|T6"D can also be 
interpreted in terms of a logarithmic spiral. Specifically, we consider the logarithmic spiral 
function m{9) defined by the equation 

m[9) = m\e = e y 6 ' ' . (17) 

The polar coordinates of the spiral are (m, 9) and m\ and a are constants. The mass 
spectrum mi, m 2 , m^, ... corresponds to the radial intercept of the spiral with the x- 
axis. The successive intercepts give the size index k. Thus m 2 = mie 27Ta and consequently 
the coefficient r = e 2na . The geometric meaning of a in terms of the logarithmic spiral can 
be obtained by considering a = d(logm)/d9 = (dm/d9)/m. The cotangent of the angle 
between any radius vector (the x-axis considered as an intercept line here) and the tangent 
to the spiral at the intercept point is then equal to a, i.e., cot f3 — a; the (3 is independent of 
the orientation of the intercept radius vector (the x-axis here) and the branch of the spiral. 
Rotation of the intercept axis (the x-axis) by angle of 9q corresponds to changes in m\ to 
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mie a9 ° and not in r, a, or (3. Secondly, it should be noted that changes in r = e 2na correspond 
to changes in a and thus in (3 since a = cot (3. As noted in Ref. ||, the logarithmic spiral is 
self-similar since m\ can be absorbed into a constant angle log m\ /a as shown in Eq. (|ITD . 

Self-similarity of the geometric sequence 1, r, r 2 , r 3 , ... generated by the coefficient r is 
also discussed in Schroeder |J. Figs. 3. 10 and 3.11 in Ref. || illustrates a triangular and a 
rectangular construction of this property of the infinite series 1 + r + r 2 + • • • = 1/(1 — r) for 
r < 1. Extending the index k to — oo < A; < oo, the series has an exact self-similar structure. 
Thus the power moments of this as considered in Ref. || has a power law behavior in 
the bin size L. However if we restrict the set to be k = 1, 2, A, then the power moment 
would have an A dependence. 

In a Malthus map the g'th order power moment of the mass distribution probability 
Pk — m-h/A with bin size L || is given by 

[A/L] 

P q (L,r)= £ 





g 


[A/L] 




E^ 




E 


r k rrii/A 


_fcej 




j=i 





'1 — r 



L\-\ 1 



„Aq 



• Lq 



(18) 



For r ~ 1, it is easy to show that P q (L,r) = (L/A) q 1 . Fig.2 shows P q (L,r) with r 
Aj (A + x — 1). Fig.3 shows the power moments of the mean mass distribution 
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k < rik > 



[A/L] 




g 


[A/L] - 


P q (L,x)= £ 


Y.<Vk> 




E 




_keJ 




j=i 



E 

fceJ 



.4 



(19) 



in x model of Ref. || compared to P q (L,r) of Eq.(^). Here, the bar over P q signifies that 
the power of the mean is being taken. Below we will also evaluate the mean of the power. 
The slope in Fig.3 corresponds to the generalized dimension D q defined in Ref. [§], 

1 log\P q (L,x) 



D q (L) 



(20) 



(g-1) log(L/A) 

This figure shows that D q is about 1 for x > 1 and much smaller than one with q dependence 
for small x. We can rewrite Eq.(|l8"D as 



P q (L,r) = P q (L,\) 



sinh(AL/2) 
sinh(AA/2) 



smh(q\A/2) 
sinh(gAL/2) 



(21) 
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where r = e A . This rapidity type quantity P q {L, A) of XL, which is normalized to 1 at L = A, 
can be used as a new variable in a study of the power moments. 

In the x model, for the bin size L = 1 and for x < 1, using Stirling approximation for 
large A » x, 



T(A + x) « (A + x - l) x - l T(A + 1), 



the power moment P q (L, x) of the mean mass distribution becomes 



(22) 



P q (l,x) 



A 

E 

k=l I 

X 



k < nk > 
~A 



A 

E 

k=l I 



x A\ T(A + x-k) 
~A(A — k)\ T{A + x) 



A 1 



[(T(x)y + (((1 - x)q)\ 



(23) 



where ((a) is a zeta function. For (1 — x)q = 1, the (((1 — x)q) should be replaced by 
\og(A — 1) + 0.57722. Since P q {A, x) = 1, then assuming the power moment follow a power 
law, we can write P q (L,x) = (L/A) aq with a q = — log[P g (l, x)}/ log(A) at low x. Then 

glog[r(l + x)] + log[l + (T(x))-%((1 - x)q)} 



On 



qx 



Usually, a q is written as 



log A 

D q (q — 1) so that the dimension D q is 



(24) 



D n 



qx 



glog[r(l + x)} + log[l + {T{x))-%((1 - x)q)} 



(25) 



q — 1 (q — 1) log A 

Fig. 4 show the power moments P q (L,x) for x = 0.5. This figure also shows a fit with an 
a q {L) = a q + b q /L form (thin solid lines) and a fit with the Stirling approximation of Eq.(f23|) 
for q = 2 (thin dash-dotted line). The D q (L) of Eq.fl2!]|) are shown in Fig. 4(b) for each q. 

Moreover, for x = m with m an integer, the mean fractional occupancy in bin J of 
length L is given by J2k < n k > I A = Pj(L,x = m), where the sum S is over k from 
k = (J — 1)L + 1 to JL. This mean fractional occupancy can be obtained from 

(A-l)\m\ 



Pj(L, x = m) 



Fj(L,m) 



(A + m- 1)! 

The Fj(L,m) = Hj_i(L,m) — Hj(L,m) and, using Eq.(48) of Ref. 

(m + Aj(L)-l)\ 



(26) 



Hj(L,m) 



(Aj(L)-l)\m\ " 
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(27) 



The Aj(L) = A - JL. At x = m = 1, the Pj{L, 1) = L/A and at x = m 
Pj(L, 2) = (2L/A) - (2 J - 1)(L/^) 2 - Th e large A Stirling limit of Pj(L, x) is 

' (J-l)LY 



Pj(L,x) 



.4 



1 



JL\ 
~AJ 



2, the 



(28) 



With this limit, the power moment P q (L,x = m) can be rewritten in terms of a simple 
convolution, namely 



P q (L, x = m) 



A 



E (Cj(m)r 

7=1 



(29) 



with n = A/L and Cj(m) the discrete convolution 

Cj(m) 



m— 1 
k=0 



(30) 



Here 6j = 1 - {J L/A). For m = 1, Cj(l) = 1 and for m = 2, Cj(2) = (n + 1 - 2J)/n. The 
geometric features in P q (L,x) arise when the factor 1 — JL/A ~ exp(— JL/A). Specifically, 

exp(-(J - l)Lx/A) - exp(-JLx/A) 



Pj(L,x) 



N(x) 



(31) 



with iV(x) = 1 — e x to guarantee that Pj is normalized to unity when summed over J. 
Then, 



P 9 (L,x) = E ^JOM) 
j=i 



sinh(ff 



sinh(g|) 
sinh(gj|; 



(32) 



sinh(|) 

which has the same structure of the geometric model of Eq . ; A <-> x/A. 

Since the Malthus model can describe only the mean behavior of the mass distribution, 
we can consider only the power of the mean distribution, Eq. (|T8|) . However for the x model, 
as we have considered in Ref. 0, we can also study the ensemble average of power moments 
in each event. For a small "integer" x, we can evaluate the power moment, 



/[A/L] 

P q (L,x) = ( £ 

\ J=l 







[A/L] I 




E^ 


') 




E kn k /A 


fee J 




J=l \ 





(33) 



of order 2 analytically. Using the relations developed in Refs. 
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[A/L] I \ V\ [A/L] 

P 2 (L,x) = £ ( *£kn k /A ) = £ EE^ < n ^k > /A 2 
j=i \ [fcg j J / j=i jcjkeJ 

_J^k 2 <n k > ^ k 2 < n k (n k - 1) > [ ^ L] jk < rijUk > 

fc=i fc=i J=i jeJk^jeJ 



(A + x) [A ^ L] ^ ^ x 2 A\ TjA-j-k + x) 
A(x+l) + hhjA 2 {A-3-k)\ 



A 2 
2 A\ 



A 2 



T(A + x) 



(34) 



Thus we have 

ACM) 

P 2 (L,2) 



1 

2 
1 
3 

W3) = ^ 



^4 + 1) L 
(A + 2) 



+ 3- 



a ' ~{a+i; 

(A + 3) 3(2A 2 + 3A-1) 



L — 



(35) 



A ' A(A + l)(A + 2) A(A + l)(A + 2) 
Eq.fl3"5|) is valid for bin lengths L which divide A into an "integer" number of parts (except 
L — A, thus L < A/2). These results show a linear dependence of the second order (q = 2) 
power moments in L rather than a pure power law or exponential behavior (see Fig.5). When 
the L independent part of P-iiL, m) is subtracted out, the resulting quantity has a power law 
behavior in L with an exponent of 1. When the constant part is accepted, these P2(L,m), 
m = 1, 2, 3, can be approximated as an exponential. In Ref. |§, we have shown that the 
power moments follow approximately the power law behavior at large x and exponential 
behavior at a small value of x. The saturation at large bin size L for large x comes from the 
high power dependences in L, such as the L 3 term in P2{L, 3). 
On the other hand, for a bin size of 1, 



Pz(l,x) 



A 

E 

k=l 



A 



A 

E 

k=l 



2 A 

< n k > + Y, 



< n k (n k - 1) > 



(A + x) + M 



A(l + x) 



k=l 



' X ' 
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A_ 





A\ 



k=l 

T(A + x - 2k) 



(36) 



(A-2k)\ T(A + x) 

Since < m k >= k < n k > is small for large k at a large value of x, then using Eq. (|22|) , the 
P2(l,x) can be approximated as 

A 



2{,) A(l + x) + 



x 



A + x -1 
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1 



A+x—l 



M^Y 



For large A » x with Stirling approximation Eq. (p2|) , 

(A + x) x x 2 T(x) 

For large A with x « A, the leading order L/A dependence in P2(L,x) is approximately 
(c.f. Eq.©) 

^ = ^71 + 1^1)+ 4' (37) 

where 

C 2 - £ + (38) 
2 A x y J 

Writing P q (L,x) oc e qL ^ q as in Ref. 0, the correlation length ( q for q = 2 is 

_ 2(A + x) 

C2 " (x + i)^' (39) 

This formula gives value for £2 which approximately agree with the results of Table 2 in 
Ref. H for small values of x. 

V. FOURIER TRANSFORMS OF MASS DISTRIBUTION AND 

THE 1/f 3 POWER LAW 

In this section we propose that Fourier spectra of fragmentation distributions be studied 
and we illustrate the procedure with an example. The importance of Fourier or power 
spectrum is the possible appearance of the ubiquitous 1/f noise which appears in many 
areas of physics. Various types of noise can be characterized by writing the power spectrum 
associated with the noise as 1// /3 , with /3 — for white, (3 = 1 for pink, (3 = 2 for brown, 
and j3 > 2 for black noise. One mechanism for generating 1/f noise is through intermittency 
generated by the logistic map near its tangent bifurcation . The possible role of the logistic 
map in nuclear fragmentation was noted in Section III. The Malthus model considered in 
this section will generate a 1/f power spectrum. This example is used for the purpose 
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of illustrating a powerful tool used in studying properties of discrete dynamical systems 
through Fourier analysis. 

The dynamical system considered here will be the set of mass points m^ for k = 1, 2, 
A or its associated probability pk = m^j A which are generated by an iterative map given by 
Eq.([Sp. The Fourier transform of m^ gives the frequency spectrum of the mass distribution; 



M f = Pk exp 



fc=i 



-i — —k 



A 



— m k exp 

k=l 



For the mass distribution of Eq.@, we have 

vn\jr " 



-i — —k 
A 



(40) 
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A 



i2irf/A 



k=l 

A J l 



-i2nf/A 



re 



-i2wf/A ' 



(41) 



Notice here ( e - i2nf/A y 



-i2wf 



1 for any integer /. An important quantity associated 



with A4f is the intensity or power spectrum; 



ml 
~A? 



1+r 



2A 



2r A cos(27r/) 



1 + r 2 _ 2r cos(2vr//A) 



A 2 V 1 



1 + (1 - cos(2tt/)) 



(42) 



1 + (i5f I 1 - cos (¥' 

The intensity If peaks at half integer / and has a valley at each integer /. For a Fourier 
transform in general, Ai f and If are periodic in / with the period of A. Furthermore, If has 
a reflection symmetry with respect to A/2, i.e., Ia-s — If- Here we consider only integer /, 
i.e., the valleys of Eq.(42) for < / < A/2. At integer /, Eq.(^) reduces to 



\M 



f\ 



mi 
A 2 



-a\ 2 r 



1 — r 



1 + 



2r 



1 — cos 



2tt/\ 



-i -l 



(43) 



(1 -r) 2 V" A 

For an infinitely large A, except for r = 1, If oc /° for finite / which is a power law in / 
with the exponent of zero. This result represent the self-similar property of the Malthus map 
discussed in the previous section. When r = 1, which is the same as the x model at x — 1, 
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Eq.(f43|) shows that // = except for values of / which are integer multiple of A. For f = A, 
I a = tu\ = 1 since m k = 1. For r very small or large compared to one which corresponds 
to a large or small value of x in the x model, the spectrum If becomes independent of the 
frequency / and thus has a white power spectrum ||; if e = r — 1 > 2irf/A, then If of 
Eq. ( f4"3"D exhibits a white spectrum. On the other hand, for r near unity, the power spectrum 
for small / compared to A can be approximated as If oc f~ 2 . The conditions for a l// 2 
power spectrum are, for e = r — 1 << 1, 

'2nf 12 



e 2 = ( r _i)2 <<; 



« 1. (44) 



.4 

In Fig.l of Section II, the choice of r ~ 1.02 fit the proton-nucleus data for A « 200. The 
eA ~ 4 and the data exhibits l// 2 behavior for 1 < / < 30 (see Fig. 9(a) in the following Sec- 
tion VI). The steepest slope occurs at / = A/ A and 3 A/4 with ± (^) 2 (^) 2 { i JT^) 2 ^f- 
As a final illustration, we consider the power spectrum associated with the result of 
Eq.([L]) at the point x=2. Substituting m k = k < n k > into Eq.(f4"0"D we obtain 

If = (A+l) 2 sin 2 (vr//A) (45) 
for / = 1, 2, . . . A — 1. With low frequencies such that irf/A << 1, If ~ l/V 2 / 2 which is 
a 1/f behavior with (3 = 2. A more detailed analysis of l//' 9 power spectrum in mass 



distribution will be given in the next section using the Au emulsion data of Ref. [10 



VI. AU EMULSION DATA 



In this section we will discuss the Au emulsion data of Ref. ||10|| . The Au emulsion data 
is an event by event analysis of the charge distribution n z , with n z the number of charged z 
fragments produced in the collision in each event. By contrast, the data considered in Fig.l 
is an ensemble average distribution. Each event conserves the total charge Z = 79 event by 
event in the Au data. 

The ensemble averaged element distribution < n z > is shown in Fig. 6 (a) and the mean 
charge distribution p z =< m z > jZ is shown in Fig. 6(b). These show the strong even-odd 



oscillation in the data (filled circles connected with a thin solid line). Here we mention that 
the charge distribution < m z > may be better than the element distribution < n z > in 
distinguishing between different model. Unlike the results of Fig.l, the Au data cannot be 
fit with a single value of x or r. The reason for this is the initial fast fall off with z in n z 
or m z . Fig.6(b) show that two mixed x's in the x model (solid line) or three mixed r's in 
the Malthus model (dash-dotted line) give a good fit to the data; both exhibit the fast drop 
of < m z > at small z followed by a flat middle range and then a rise at the large z. Also 
shown in Fig. 7(a) is the multiplicity distribution P(M) where the multiplicity M = J2z n z, 
event by event. This figure shows that the emulsion data does not corresponds to a single x 
nor a single r. 

In Refs. the importance of the cumulative mass distribution 

z 

M(z) = M Z = Y, rrij (46) 

3=1 

was also stressed. The M z has a staircase behavior and the sequence M 1? M 2 , M z , 
M z can be considered as an ordered sequence similar to the sequence of "energy" levels || . 
If we use the mean distribution < m z > in Eq.([5J), then, at x — 1 or r = 1, levels are 
uniformly spaced with a spacing s = 1. The distribution of the experimental "spacings" 
is then determined by the experimental values of = kn^ event by event which may 
then be compared with a Poisson or a Wigner distributions. Fig. 7(b) shows the m z spacing 
distribution in the Au emulsion data (solid circle connected by a thin line). Here m z is 
the spacing in each event and the probability is the ensemble averaged one. The data looks 
more like a Wigner distribution (m 2 e~ mz//5 /25) than a Poisson spacing (e _mz//10 /10). In level 
spacings, the Wigner distribution is discussed in terms of a chaotic structure in the spacing 
distribution. It should also be noted that the staircase feature || of the cumulative mass 
distribution of the data (thin solid line in Fig. 8 (a)) can be described by two x's (thick solid 
line) or by three r's (dash-dotted line). These three curves are quite similar as can be seen 
in the figure. 

In Fig. 8(b) we show a cobweb type plot (m z vs m 2+1 ) by solid circles. Also shown in 
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the figure are the pairs of (m z , m 2+2 ). The origin of the scatter of the points is the even- 
odd oscillating behavior which also has a varying amplitude in the charge distribution of 
Fig.6(b). If the data follows the x model or the Malthus model, then the points would lie 
on a single line. The scattered points indicate that even the mean distribution of the data 
exhibits some random structure. 

Finally, we show the Fourier spectrum of the data in Fig.9. Fig.9(a) shows the power 
spectrum of the data compared with the power spectrum of the x model and the Malthus 
model. In the low frequency region, the x model with x = 0.3 (dashed line) exhibits a 
behavior with (5 ~ 0.56, the Malthus model with a single r (dotted and dash-dot-dot-dotted 
lines) has (3 xs 1.9, and the Malthus model with three mixed r's (dash-dotted line) follows 
a power law with (3 xs 0.68. The data (solid circle connected by a thin line) shows a more 
complicated structure than the x model or the Malthus model. Also shown by open circles 
connected with a thin dashed line is the Fourier spectrum of the data in which the high 
peak at z = 1 and 2 in the m z distribution is reduced by putting these values at the value 
of m 3 . This spectrum then follows the single Malthus model at the low frequency part. In 
Fig.9 (b), we show the / dependence of the fluctuations of the data from the mean behavior 
given by the x-model. The thin dashed line is for x =3, the thick solid line is computed 
with 2 a;'s (66% x = 0.5 and 34% x = 50), and the thin solid line is the Malthus model with 
three r's. This figure does not exhibits a simple 1/ f 13 behavior. 

VII. CONCLUSION 

In conclusion, this paper illustrates a simple approach to the distribution of mass pro- 
duced in a heavy ion collision. This distribution can be considered in terms of iterative one 
dimensional maps similar to those used in population dynamics. In this picture, the orbit 
of a monomer or nucleon in its approach to an attractor of the dynamics generates the mass 
distribution. Simple Malthusian dynamics seem to account for some features of the data. 
A non-linear model is also discussed to see if the data has any evidence for non-linearities 
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and limiting behavior in otherwise geometric growth. When the Verhulst model is used to 
generate the mass distribution m^, the resulting distribution of cluster sizes rik = m^jk is 
shown to have an interesting bifurcation point whose signatures are discussed. New ways 
of representing data which go beyond any specific model are also suggested. These ways 
include plotting pairs (m fc , m k+1 ) which appear in cobwebbing, plotting mass yields, cumu- 
lative mass yields, and inverse mass yields versus cluster size k besides the traditional yields, 
n>k = rrik/k versus k. More complex maps can also be investigated than those discussed here 
as well as fluctuations such as in Ref. ||. 

Using the simple geometric model developed in this paper, questions related to self- 
similarity and scaling behavior are discussed. Simple scaling relations for the power moments 
of the mass distribution are developed based on self-similar geometric properties of the mass 
yields. The correlation lengths, exponents and dimensions associated with these scaling 
relations are discussed. The power spectrum associated with Fourier transforms of the mass 
spectrum is also developed. Under certain conditions a 1/ f 13 spectrum is found. 
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FIGURES 

FIG. 1. Fragment mass distributions in the radiochemical measurements from the decay of a 
target residue produced in a collision of high energy proton on a nucleus. The filled circles are the 



data summarized in Ref. [12] and the histogram are the thermodynamic statistical fits of Ref. [12]. 
The solid curve is obtained using the x- model of Eq.(|l|) with x = 0.3 and the dashed line is the 
orbit under the Malthus linear iterative map of Eq. @ . The coefficients r for the Malthus map are 
r = 1.065 in (a), 1.027 in (b), 1.024 in (c), and 1.041 in (d). These values are the geometric mean 
of the corresponding rA_(k,x) with x = 0.3 in all four cases. The dash-dotted line in (c) is a linear 
map with the recursive coefficient = 1 + (k — 50)/3000. The dash-dotted line in (d) is an orbit 
in the non-linear Verhulst map Eq.© with C = 0.062, m s = 170, and mi = 1. 

FIG. 2. The power moments P q {L,r) of the mass distribution in a Malthusian model. The 
mapping coefficients r for each x are r = A/ (A + x — 1) with A = 79. The solid, dashed, dotted, 
and dash-dotted lines are for q = 2, 3, 4, and 5 respectively. 

FIG. 3. The power moments P q {L,x) of the mean distribution in the x model vs the power 
moments P q (L,r) of the Malthus model of Fig. 2 with the corresponding x values. The lines are 
same as in Fig. 2. 

FIG. 4. Power moments P q (L,x) of the mean distribution of mass in the x model for x = 0.5 



and A = 79 (a) and their corresponding slopes, a q (L) = (q — l)D q (L) using Eq.(^0|) (b). The solid, 
dashed, dotted, and dash-dotted thick lines are for q = 2, 3, 4, and 5 respectively. The thin solid 
lines in (a), which are close to the thick lines, are fits using a form P q (L) = (L/A)( aq+bq / L )( q ~ 1 ^ and 
the thin dashed lines are the simple power law fit of P q {L) = (L / A) Dq( - q ~ 1 \ The thin dash-dotted 
line, which is close to thin dashed line, is based on Stirling's approximation for q = 2, Eq.( 
Here, a q = D q - b q with b q = 0.04 and D x = 0.942, D 2 = 0.849, D 3 = 0.756, D 4 = 0.695, and 
D 5 = 0.656. 
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FIG. 5. Power moments P q (L,x) which are the ensemble average of the power of the mass 
distribution in x model for A = 79. The solid and dashed lines are for q = 2 and 3 respectively. 
The dotted lines for x = 1.0 and 3.0 are the exact results of Eq.(|35"|) and the dotted lines for x = 0.5 
and 5.0 are the approximation of Eq.(|37|). 

FIG. 6. The cluster distribution (a) and the charge distribution (b) of Au emulsion data 
(Z = 79) of Ref. [10]. The solid circle connected by a thin line is the data. The dashed line 



is the x model with x = 0.3 and the solid line is the x model with a mixture of 0.66 of x = 0.5 and 
0.34 of x = 50. The dotted line is the Malthus model with r = Z/(Z + x — 1) with x = 0.3 and the 
dash-dot-dot-dotted line is for r given by geometrical mean of x model; the corresponding r values 
are 1.0089 and 1.0546 respectively. The dash-dotted line is obtained from a mixture of three r's: 
1.1210 (29%), 0.9938 (40%), and 0.3492 (31%) respectively. 

FIG. 7. The multiplicity (a) and the "level spacing" (b) distributions in the ensemble of 415 
Au data events (the solid circle connected by a thin line). In (a), the dashed line is the x model 
with x = 0.5, the solid line is the mixture of 0.66 of x = 0.5 and 0.34 of x = 50, and the thin solid 
line is the mixture of 0.72 of x = 0.5, 0.15 of x = 5, and 0.13 of x = 79. The dashed and solid lines 
in (b) are Poissonian (e _mfc//10 /10) and Wigner (mj.e~ mfe//5 /25) spacing distribution respectively. 

FIG. 8. Staircase (a) and cobweb (b) plots. In (a), the curves are the same as in Fig. 6 except 
the thin solid line is for the data which is very close to the dash-dotted line and to the thick solid 
line. In (b), the solid circles are the pairs of (m z , m z +\) and the closes are the pairs of (m z , m z +2)- 
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FIG. 9. The Fourier transform of the charge distribution p z (a) and the noise spectrum (b) for 
the Au data. The curves in (a) are the same as in Fig.6. The open circles connected by a thin 
dashed line, are the Au emulsion data except that the charge distribution at z = 1 and 2 are each 
taken the same as the z = 3 value. This is done to reduce the effect from the large probability of 
charge number 1 and 2. In (b), the thin dashed and the thick solid lines are the noise spectrum 
of the Au data with respect to the x model with a single x (x = 0.3) and for two mixed x's (0.66 
of x = 0.5 and 0.34 of x = 50) respectively. The thin solid line is the noise spectrum with respect 
to the Malthus model with three r's as in Fig.6. Also shown by the solid circles connected with a 
thin line is the power spectrum of the data which is shown in (a). 
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